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We study the dependence of the free energy on the CP violating angle 9, in four-dimensional SU(iV) gauge 
theories with N > 3, and in the large- TV limit. Using the Wilson lattice formulation for numerical simulations, 
we compute the first few terms of the expansion of the ground-state energy F(9) around 9 — 0, F(9) — F(0) = 
A 2 e 2 (l + b 2 9 2 + ...). Our results support Witten's conjecture: F{9) - F(0) = AO 2 + 0(1/N) for 9 < ir. We 
verify that the topological susceptibility has a nonzero large-iV limit \oo = 2.4 with corrections of 0(1/N 2 ), 
in substantial agreement with the Witten-Veneziano formula which relates Xoa to the Tj mass. Furthermore, 
higher order terms in 9 are suppressed; in particular, the 0(8 4 ) term b 2 (related to the if — rj elastic scattering 
amplitude) turns out to be quite small: b 2 = —0.023(7) for N — 3, and its absolute value decreases with increasing 
N, consistently with the expectation b 2 = 0(1/N 2 ). 



1. INTRODUCTION 

Interest in SU (N) gauge theories has recently 
been revived by studies in the context of M- 
theory, AdS/CFT correspondence, and N = 2 
SUSY broken to JV = 1 . A number of predictions 
have been put forth, and we have set out to check 
these predictions by numerical simulation. 

One aspect of these studies regards the spec- 
trum of confining strings in four-dimensional 
SU(A^) gauge theories. A considerable amount 
of numerical work has been done on this sub- 
ject (see and references therein). In par- 
ticular, in Rcf. Q wc have presented computa- 
tions of the string tensions , related to sources 
with Zn charge k, for SU(4) and SU(6) gauge 
theories. Our results are consistent with the sine 
formula a^ju\ — sin(fc7r/iV)/ sin(7r/7V), with an 
accuracy of approximately 1% [2%] for SU(4) 
[SU(6)]; this formula emerges in various realiza- 
tions of supcrsymmetric SU(-/V) gauge theories, 
and in the context of M-theory. On the other 
hand, our results show deviations from Casimir 
scaling a^ja = k(N — k)/(N — 1), as might be 
expected from the short distance behaviour of the 
potential, or from a strong-coupling limit. 

In the present work, we study the 9 depen- 
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dencc of SU(iV) gauge theories, for N > 3 and 
in the large- N limit, by numerical simulation; for 
a longer write-up of our work, sec ||. The angle 
9 appears in the Euclidean Lagrangian as: 

C e = \F; v {x)F- v {x)-z9q{x) (1) 

where q(x) is the topological charge density 

q{x) = 64^ W^(z)^U*)- (2) 
The partition function reads: 

Z{9)= J[dA]exp(~ Jd 4 xC e ) = exp[-VF(9)} (3) 

where F(6) is the free (ground state) energy. Wit- 
ten has argued that in the large- N limit F(9) 
is a multibranched function of the type 

F(9) = iV 2 min fc H[{9 + 2nk) /N] (4) 

The conjecture was refined [pi and simplified: 

AF{9) = Amm k {9 + 2nk) 2 + O (1/N) . (5) 

(AF(9) = F{9) - F(0)). In particular, for 9 < tt, 

AF(9) = A9 2 + 0(l/N). (6) 

This conjecture has been supported by 
gauge/string duality arguments Q. 
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Monte Carlo studies of the ^-dependence are 
made very difficult by the complex nature of the 
9 term. In fact 6^0 cannot be directly simu- 
lated on the lattice. Here we restrict ourselves to 
relatively small 8 values, where one may expand: 



AF(6) = A 2 9 2 + A±6 4 



(7) 



A 2 gives the topological susceptibility, i.e. 
^2=x/2, x = <Q 2 W^ (8) 

where Q is the topological charge Q = J d 4 x q(x). 
Higher order coefficients in Eq. can be related 
to higher moments of the probability distribution 
p(Q) of the topological charge. For instance, 



A 4 = 



2W 



[<Q 4 ) - 3(Q 2 



9=0 



(9) 



A convenient parameterization of the free energy 
using dimensionless scaling coefficients is: 



1 



<j- A AF(6) = -CP (I 



b 2 e 2 + b 4 e i + ...) (10) 



(a: string tension, C = x/c 2 , b 2 = A4/A2). 

Witten's conjecture implies that C has a finite 
nonzero large- TV limit Coo and b 2 i = 0(1/N 2 ). In 
particular one expects that b 2 is small, and that 
it should rapidly decrease with increasing N. A 
nonzero value of Xoo is essential to resolve the 
U(1)a problem, and is related to the 77' mass 0: 



Xoo = 2A 



ANt 



Q(l/N). 



(11) 



2. CALCULATION AND RESULTS 



We present results for four-dimensional SU(iV) 
gauge theories with TV = 3,4,6. In a nutshell, 
we: Determine the ratio C = x/ " 2 ; verify the 
expected behavior C ~ A + B/N 2 , determining 
the large- N limit by a fit; obtain a rather accurate 
estimate of Coo; estimate the values of b 2 , which 
are found to be very small and rapidly decreasing 
with N, in support of Witten's conjecture. 

Simulation Details: • We use the standard 
Wilson gauge action and the Cabibbo-Marinari 
updating algorithm (one SU(N) update — ► N(N— 
l)/2 SU(2) subgroup updates). • Microcanonical 



over-relaxation and heat-bath updates are used in 
a 4:1 ratio. • Values for (3 are chosen in the weak- 
coupling, (beyond the 1st order phase transition 
(N = 6), beyond crossover (peak of specific heat 
.ZV = 3, 4)). • Topological quantities are com- 
puted via cooling. We measure Q typically every 
100 sweeps. • We have employed different lattice 
sizes, so that finite size effects are under control. 
Tables of our numerical results can be found in 
Ref. [i. 




FIGURE 1 . Autocorrelation time tq versus £o- = 
a -1 ' 2 . Dotted lines show the linear fits. 

A severe form of critical slowing down is seen in 
the measurement of Q, at large values of N. Our 
results (Fig. 1) suggest an exponential behavior: 

In tq w N (e N & + c N ) (12) 

= a -1 / 2 , and ejv, c/v are constants with only 
weak A-dependence) . Thus, In tq oc N (at fixed 
£<j). This phenomenon is possibly due to the free- 
energy barriers present. There are no a priori ar- 
guments in favor of this behavior; it simply arises 
naturally from fits, unlike a power law. No simi- 
lar effect is observed in plaquette or Polyakov line 
correlations; this suggests a decoupling between 
topological modes and nontopological ones, such 
as those determining confining properties. 

In Fig. 2, the ratio C is plotted versus cr, to 
evidentiate possible scaling corrections (expected: 
0(a 2 ), logs). Corrections are indeed observed for 
N = 3 and N = 4, and we must allow for them 
in a fit, for a continuum extrapolation. Assuming 
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a linear behavior: C — C con t + ua, to take into 
account the leading 0(a 2 ) correction, we obtain: 

C cont = 0.0282(12), u = 0.077(24) (N=3) (13) 
Ccont = 0.0257(10), u = 0.049(15) (iV=4). (14) 

For N = 6, there is no evidence of scaling correc- 
tions; we find: C cont = 0.0236(10) for N = 6. 
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FIGURE 2 . The scaling ratio C — x/ a versus a. 
Dotted lines show the linear fits. 

In Fig. 3, our results for C con t are plotted ver- 
sus 1/N 2 (the expected order of 1/7V corrections). 
They are clearly consistent with a behavior of the 
type Ccont (N) = Coo + B/N 2 . Assuming such a 
behavior, a fit to N = 3, 4, 6 results gives 

Coo = 0.0221(13), 5 = 0.055(19), (15) 

This is a rather accurate estimate of large- V limit 
of C, and it is clearly nonzero. 

To compare with literature, we set a standard 
value: y/a = 440 MeV, for N > 3. We obtain: 
X H 4 = 170(3) MeV, and x 1/4 = 180(2) MeV for 
N = 3. From Eq. JTl]), with actual values of 
f„, m v i, and Nf — 3, one has 191 MeV. Using 
Veneziano's improvement: m 2 , — ► m 2 , + m 2 — 
2m 2 K7 one obtains 180 MeV. 

We now turn to 62: The cancellations in its def- 
inition necessitate very high statistics. We have 
obtained sufficiently accurate results for N = 3, 4, 
showing scaling. A weighted average, with con- 
servative errors, yields: 

b 2 = - 0.023(7)(iV=3), 6 2 = - 0.013(7)(JV=4).(16) 



Reducing statistical error on b 2 in the N = 6 
case is extremely difficult. We can only give a 
rough estimate b 2 = -0.01(2) (7 = 0.348) in this 
case, showing that b 2 is very small. 
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FIGURE 3. The scaling ratio C = x/°" 2 versus 
1/N 2 . The line is the linear fit to A + B/N 2 . 

An alternative to cooling is based on the in- 
dex theorem, using, e.g., the overlap formulation 
M. The rigorous properties of such an approach 
would make up for a useful consistency check. 

In conclusion, we stress that further investi- 
gation of other observables in SU(N) theories 
is called for, to test predictions from: SUSY 
SU(N), M-theory and gauge/string duality. 
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